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We study the strange quark mass dependence of the tricritical point of the [7(3)1, x f/(3)ij linear 
sigma model in the chiral limit. Assuming that the tricritical point is at a large strange mass value, 
the strange sector as well as the r) — ao sector decouples from the light degrees of freedom which 
determines the thermodynamics. By tracing this decoupling we arrive from the original U{?>)l x 
U{?i)r symmetric model, going through the U(2)l x U{2)ii symmetric one, at the SU{2)l x SU{2)ii 
linear sigma model. One-loop level beta functions for the running of the parameters in each of these 
models and tree-level matching of the coupling of these models performed at intermediate scales are 
used to determine the influence of the heavy sector on the parameters of the SU{2)l x 5[/(2)_h linear 
sigma model. By investigating the thermodynamics of this latter model we identified the tricritical 
surface of the U{3)l x U{3)r linear sigma model in the chiral limit. To apply the results for QCD 
we used different scenarios for the rris and fiq dependence of the effective model parameters, then 
the fj,]i (ms) function can be determined. Depending on the details, a curve bending upwards or 
downwards near /i, = can be obtained, while with explicit chemical potential dependence of the 
parameters the direction of the curve can change with nis, too. 



^ . I. INTRODUCTION 



The phase diagram of QCD is a much-studied phenom.enon, but still its characteristics, at finite baryon chemical 
potential in particular, are far from being settled [1-3]. While at zero chemical potential all the Monte Carlo (MC) 
and effective model studies tend to support a common picture, at nonzero chemical potential the MC results are 
inconclusive. 



> 

QQ : At zero chemical potential the widely accepted phase diagram in the mud-iTT-s plane exhibits both at small and at 

fS| large quark masses regions of first order phase transition each bounded by a line of second order critical endpoints 

(CEP). In between these regions, the transition is of analytic crossover type. If we introduce a nonzero quark baryon 
(^ •■ chemical potential fiq, the second order CEP lines extend to a critical surface. If the critical surface lying closer to 
f^ [ the origin of the mass-plane bends upwards, that is to larger quark masses, then there is a possibility that a crossover 
transition becomes at larger chemical potential a second, and subsequently a first order phase transition. Direct MC 
simulations [4, 5] and also the estimates in [6, 7] give a finite n^^^ value, and similar conclusions can be drawn 
from other lattice techniques, too [8-10]. However, the second order surface seems to bend downwards, to smaller 
^^ quark masses, according to a lattice study of the curvature performed in [11-13] using imaginary chemical potential. 

Although all these lattice results do not necessarily contradict each other, they could imply a scenario in which the 
critical surface has a non-trivial shape. Some numerical evidence for such a possibility is given in [14]. 

Beyond direct simulations one can approach the study of the phase structure of QCD through effective theories, 
see e.g. [3, 15] for reviews. At finite chemical potential there are several results on the chiral phase transition for two 
flavors in Nambu-Jona-Lasinio (NJL) models, as well as in linear sigma models [16-19]. Using the SU{2)l x SU{2)fi 
linear sigma model an interesting phase structure with two CEPs in the fiq — T plane was reported for low values of 
the pion mass in [19]. Considering three flavors, one can study in these models the properties of the chiral critical 
surface [20-25]. In [20-22] the authors used the U{3)l x C/(3)_r chiral sigma model near the physical point, and found 
that in the available parameter space the critical surface bends upwards, supporting the direct MC result. In [23-25], 
in the framework of the extended NJL model the authors found a down-bending surface for small chemical potential 
which eventually turns back at higher values of /ig. This behavior would conciliate the two MC scenarios within a 
single critical surface if the turning happens at positive values of mud, fn^ and finite values of iiq. Other possible 
behaviors of the critical surface were discussed in [26, 27] using the Gibbs' phase rule for phase coexistence. 
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In the chiral limit {rriud = 0), there are two weh- known hmits, the fiq = and the infinite strange quark mass 
{mg — oo) hmits. At /i^ = we have a first order phase transition region for smah rUs, a second order transition region 
at large nis and a tricritical point (TCP) separating them, with a characteristic value mj'~^^ . At nis — oo and small 
chemical potential we have second order, for large fiq a first order phase transition, with again a TCP in between. 
The line of TCPs is the intercept of the critical surface and the mg — fJ.q plane. According to the above scenarios the 
two TCPs at the /iq = and mg ~ oo can be connected by a single line with definite curvature, a back-bending line, 
or they may belong to two distinct TCP lines. In the latter case there must be two disjunct critical surfaces in the 
TTiud — TTis — fJ'q spacc, of which these two TCP lines are just the endpoints in the mg — iiq plane [26, 27]. 

In this paper we attempt to describe the behavior of the tricritical line in the chiral limit of the U(3)l x U{3)r 
sigma model [20, 28, 29], assuming in addition that the mass of the constituent strange quark and the anomaly scale 
are much larger than the critical temperature. The study of the chiral limit has some advantages, since one can work 
with much less parameters than in a generic situation. A disadvantage though, is that in this case there are no direct 
measurements which could connect the effective model parameters with the QCD (although, strictly speaking, there 
are no such measurements anywhere, apart from the physical point, especially not for a possible chemical potential 
dependence of the parameters). The goal of this study is to explore the parameter dependence of the TCP line. 

The assumption that the constituent strange mass is much larger than Tc is based on the observation that even at 
the physical point of the mass plane the constituent strange quark mass is ^ 450 MeV, while Tc ^ 160 MeV, and the 
critical line in the rriud — ms plane behaves as mj^^ — rUg ~ m^^ , as one approaches the rUg axis, which predicts for 
toJ^^ and for the mass of the strange constituent quark a higher value than the corresponding mass at the physical 
point. Actually, in the lattice study of [30] it was estimated that mri^'^^ — 2.8Tc, while in [29] using the U{3)l x U{3)b, 
sigma model to J was estimated to be one order of magnitude bigger than the value of nig at the physical point. 
A similar observation can be made for the anomaly scale, which is connected to the mass of the 77' meson. In this 
physical situation we have a multi-scale system, where a simple one-loop analysis would lose the contribution of the 
heavy sector. Instead, we have to work with decoupling theory [31] which results in a hierarchy of effective models 
describing the physics at lower and lower scales: first the mg strange quark sector decouples with the corresponding 
bosonic degrees of freedom, and we obtain an effective U{2)l x t/(2)/j symmetric theory. Then the 77' sector (which 
has semi-large masses at m^ — > 00 because of the anomaly) decouples, and we are left with the SU{2)l x SU{2)r 
chiral sigma model consisting of pion and sigma mesons, as well as the u and d constituent quarks. Here we can use 
the results of [18] to determine the position of TCP for a given parameters set. The effect of the strange sector on 
the position of the TCP is only through the modified parameters of the SU{2)l x SU{2)ji chiral sigma model. To 
determine the parameters of the different effective models involved in the analysis we use the one-loop /3-function 
governed running of these parameters and matching of the corresponding n-point functions in the common validity 
range of the models. Some extra assumptions on the original parameters of the U{3)l x 1/(3)^ sigma model are 
unavoidable. 

The setup of the paper is as follows. First, we discuss the model in Sec. II. Then, we overview the generic ideas 
how the decoupling works in Sec. III. Next, we perform the decoupling in the U{3)l x 1/(3)^ model in Sec. IV. We 
study the thermodynamics in the resulting effective SU{2)l x SU(2)ji linear sigma model in Sec. V. We close with 
conclusions in Sec. VI. 



II. THE MODEL 

We first construct the starting model [20, 28], which is the U{3)l x U{3)u symmetric linear sigma model defined 
by the Lagrangian: 

Cui3) =^i$- 2gT^aa + i"f57ra)] ib + Tr(9^$ta^$) _ [/($)^ (1) 

where T" are the U{3) generators, Ta = Xa/2 for a = 1, ... 7 and T° = T^, T^ = TV [29], with 

1 /ioo\ 1 /ooo\ 

T^ == - 1 r« = ^ . (2) 

^ yoooy v2 yooiy 

The meson matrix $ = cr + ztt = T°'{(Ja + iT^a) is given in terms of the physical degrees of freedom as: 

M^.+al) a+ K+\ ^ /i(^^ + ^0) ^+ ^4 




}^,ia. -«!])«:", ^=^\ ^~ ^,iv.^ ^°) K^ 1 • (3) 



Finally, the potential for $ reads 

t/($) = M2$t$ + Ai [ Tr($^$)] V A2 Tr [(^'f^)^] - \/2C(det $ + det ^l') + Tr [i7($ + ^l')] . (4) 

Here, C governs the U{1) anomaly which breaks the symmetry to SU{3)l x SU{3)u. The last _ff-dependent term 
explicitly breaks also this symmetry. As a consequence we assume that vacuum expectation values develop for those 
scalar fields belonging to the center elements of the symmetry group. These are taken into account through the shifts: 
ax -^ Gx + X and CTj, — > (Ty + y. In this paper we are interested in the regime where the vacuum structure of the 
theory contains a heavy s-quark sector and a light ud sector, with a corresponding meson sector. For the constituent 
quark sector this mass hierarchy can be fulfilled with y 3> x, since at tree-level the s-quark mass is to^ — v25J/7 while 
fTiud = gx. We assume that the mesons containing s-quark have a mass of order mg, while the rest have a mass of 
order vriud- No splitting is assumed between u and d quark masses. 

After having made all these assumptions, the mass spectrum, and correspondingly the physics splits into a heavy 
and a light part. The light sector influences only very little the symmetry breaking pattern of the heavy sector, 
therefore we may set a: = when dealing with the determination of the heavy masses. Then we find the following 
mass relations: 

ml = M^ + \iy^-Cy, ml ^ M^ + X^y^ + Cy, mj^ ^ M^ + (Xi + X2)y^, 

ml = M^ + i{X^+X2)y\ m^ = 2gV, (5) 

where m^ is the mass for tt^ and Gx, m^ is the mass of rjx and a°, m\ applies for riy, K and k modes, m^ is the mass 
of (Tj, and finally, nig is the strange quark mass. In case of a large strange quark mass, i.e. large y values, the lightness 
of the pion-sigma sector requires that rn^ <^ mj.. This can be achieved only with fine-tuning rn^ ~ after all this is the 
manifestation of the hierarchy problem, where the high energy sector influences, through radiative corrections and 
spontaneous symmetry breaking, the light sector. To circumvent the problem, we parametrize everything with the 
light x-mass: 

ml^ml+2Cy, m% ^ ml + X2y^ , ml = ml + {2Xi + 3X2)y^. (6) 

We see that there is a third mass scale, associated to Cy. Since C has the dimension of a mass, we can relate its value 
to y or m,x, and it is a matter of the low- lying dynamics determining the details of the effective model to know which 
is the "true" ratio between them. In this work we try to play around the possible values of C to see its effect on the 
thermodynamics. 

So, after all, we have three possibly very different mass scales, mg, vria and nix. The thermodynamics of the system, 
which is related to the spontaneous breaking in the non-strange sector, must take place at the light scales, that is 
at r ^ nix- To treat a physical system with vastly different mass scales is possible only using the fact that for 
the light physics the heavy degrees of freedom decouple, and their presence can be identified through the values of 
the parameters of the Lagrangian containing the light degrees of freedom. How this decoupling works in detail, is 
summarized in the next section. 

III. DECOUPLING OF MASS SCALES 

Here we review the generic principles of decoupling described in details in [31]. Let us start with a theory with 
Lagrangian £, couphng constant set {gi}i=i...u (this includes also the masses) and field contents generically denoted 
with $ and Lp. Let us assume that $ is much heavier than (p, their masses are denoted by M and to, respectively^. 
Then we want to establish an effective theory based exclusively on the light degrees of freedom. Let us denote the 
Lagrangian of this effective theory by £, its couplings are {gi\i=i...u [u < u) and the field content is (f with mass 
rh ^ m. We emphasize that, although ip corresponds to the light degrees of freedom (p of the complete theory, they can 
differ by wave function renormalization. In general the wave function renormalization and the renormalized couplings 
of the effective theory depend on the coupling constants of the original theory: 

^^zy\g)^, m^Mg). (7) 

If we assume that the couplings of the effective model at tree level are linear combinations of the original couplings, 
then we can write 

u 

9t=^Yl ^'^9] + Affj (.9) , i€[i...u\, (8) 



For the sake of simplicity we only consider one heavy and one light degree of freedom 



where Agi{g) denotes the loop corrections which depend on aU of the original couplings. 

In order to determine these relations we use matching: we take u + 1 correlators of the light fields at some external 
momentum, and require that their physical value be the same in the original and in the effective model. If we denote 
the n-point functions of the original and effective model by 



G'"(xi,...X„) = {T(p{xi)...(p{Xn)) , G"'{xi,...Xn) = (T(^(xi) ...(p{Xn)) , 



then we require in Fourier space 



G"(fci,...fc„) = z"/2G"(fci,...fc„) 



(9) 



(10) 



for fixed fci, . . . fc„. The right hand side contains, as radiative corrections, the effect of heavy modes. 

In perturbation theory the n-point functions depend also on the renormalization scales /i and fi, respectively. In 
general, we expect to obtain all type of logarithmic corrections In" fi/E, where E is any energy scale which shows up 
in the given n-point function. In loop integrals we typically find E"^ ~ max(fcf ,mass^), where the mass can be m, M 
on the right, or m on the left hand side. In order not to have multiple mass scales in the n-point functions (which 
would lead to large logarithms) we shall choose \ki\ ~ M, then all scales are of the order of M. This means that for 
the best convergence of the perturbative series we shall also choose /x, /i ~ M. In consequence the relation (7) is in 
fact established at scale M: 



,(M) = J2 Sv9AM) + Ag,ig{M)). 



(11) 



In the original theory the couplings may be defined on a different energy scale, Mq. Then, we have to run the couplings 
according to C in order to find gi{M) which is needed above. On the other hand, (II) defines cji on the scale M . If 
we need them on a lower energy scale M', then we have to run them according to C, the effective model Lagrangian. 
If there are multiple scales to decouple, we have a series of effective models C^^'' with coupling constants {^j }i=i...M„ 
where the heaviest field has a mass M(^a) ■ Then the matching conditions described above lead to the series of equations 



9i 



(a+l) 



[Mi^,)) ^Y.^,m\M(a^) 



A5f^(ff(M(,))). 



(12) 



This defines the new couplings g^ at the mass scale Mf^^-^ which, from the point of view of the a-l-lth effective model, 
is a high energy scale. Then the running of the couplings between Mi^x — > M(„_|_i) is governed by the Lagrangian 
£(a+i). This process leads to the schematic running depicted on Fig. 1. 
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FIG. 1. Schematic plot of the running of a coupling constant. The running between Af(a) and M(a+i) is governed by C(a+i)- 



In the most simple version we use tree level matching, which means that we neglect the Ag^ terms. We use this 
approximations throughout this work. Moreover, we use the lowest order (one-loop) beta functions for running. 

We start with the U{i)L x C/(3)_r symmetric linear sigma model which is parametrized at the heavy, strange mass 
scale in order to avoid the need for running. After the decoupling of the strange sector we are left with a,U {2) lxU [2) fi 
model, starting at vris scale. The beta- functions of this model determines the running of the coupling down to ?7ia, 
when wc switch to SU{2)]^ x SU{2)fi model containing only the a and tt mesons and the u and d constituent quarks. 
We perform in this model the running of their couplings from rria down to T scale, where the finite temperature study 
is performed. 



A separate treatment is needed when we treat the light masses of the system. At one hand, due to the RG running, 
they acquire logarithmic dependence on the heavy scale. On the other hand, because of the hierarchy problem 
discussed in the previous section, the light mass squared rn^ have to be fine tuned by M^ in order to keep the light 
sector truely light. Therefore in bosonic models the logarithmic corrections from the RG running are subleading, and 
so they should be neglected. 

IV. DECOUPLING OF THE HEAVY SECTOR 

We use (1) with a background y taken into account via the shift ctj, ^> CTj, + y. The first step is to parametrize 
the heavy sector by determining the parameters of the model at scale m^. Then, as a next step, we perform the 
decoupling of the heavy (strange) sector. 

A. Parametrization of the heavy sector 

We should fix the parameters of the heavy sector by measurements, but there exist no direct mass measurements 
in this regime. Therefore we are forced to make some assumptions. Since we arc in a large quark mass regime, we can 
not use the results of the chiral perturbation theory. Instead, the heavy constituent quark model approach, where the 
mass of the heavy particles are simply the sum of the constituent quark masses, seems to be more adequate. In case 
of the light-heavy mesons this works nicely, since we can require 

niK — TTLs => \2 = 2g . (13) 

In the doubly heavy sector there is a mismatch between the scalar (ay) and pscudoscalar {rjy) meson masses, although 
the constituent quark model would give the same mass for both. The reason in this model is that there is a symmetry 
breaking effect in addition to the constituent quark masses. To treat this situation we introduce a free parameter A, 
and require that some average of the two mass squared is (2™^)^: 

Aml + {1-A)ml^i2m,sf, (14) 

which results in a relation between the two quartic coupling constants Ai and A2 

A, ^ (^ - 1) A. (15) 

If ^ = 1, i.e. when iriy — 2?Tis, Ai == g^ is the smallest value for Ai. Another plausible choice is ^ = 1/2, then 
Ai = 2A2 = 4g^. Therefore, introducing A = —2 + 3/ A we may set 

Xi=Ag'^, where Ae[l,oo]. (16) 

The parameter range A (z [1,4] will be used later in the analysis. 

For fixing the other parameters we will use the infrared (IR) sector: we will determine m^ at its own scale, this 
will give the mass unit in the study. For fixing the rris value at /i = rris scale we also use an IR observable: this will 
be that rris value where we have a TCP at zero chemical potential (but, of course, at finite temperature). Through 
the decoupling equations this value will determine the original, zero temperature strange quark mass. 

B. Tree level decoupling 

The next step is to eliminate the strange sector, and determine the parameters of the resulting effective model. The 
effective model will contain the following degrees of freedom: the upper two (m, d) components in ip, denoted by -02, 
and the upper left 2x2 submatrix in $, denoted by $2 and containing the a = a^, oq, rjx, and tt fields: 

^-ft). ^^fi!:, i^"\V — -4:^^:5?) ■ <") 



We split the Lagrangian (l)-(4) according to this characterization of degrees of freedom. The shifted Lagrangian reads 

Cu(Z) = Cu(2) + >Choavy + COnSt. . (18) 

Here, Cheavy contains the heavy part, the 'const.' term refers to the y-dependent part of the potential, while ^u(2) 
contains the terms which consist of the light fields: 

>C[/(2) =^2 [i0 - 9T^{(J2^ + ^TsTTz^)] ^2 + Tr (9^$^9^$2) - (M' + Aiy^) Tr [$t$2] 

-Ai [ Tr($J$2)] ^ - A2 Tr [($^$2)^] + Cy [dot $2 + dot $^1 . (I9) 

Tree level matching means that we simply neglect the heavy part, and go on with the Lagrangian described above. 
We can argue for this simple choice as follows. The parametrization of the heavy sector of the model could be done 
only very heuristically, in which case only the leading order effects could be taken into account. Therefore, it would 
be inconsistent to work with a detailed decoupling scheme, determining Ag ^ 0{g^) or 0{g^) corrections to the tree 
level values which are not known precisely. For this reason the leading order approach is the most consistent here. 
As a result of this approximation we shall trust only the most robust consequences of this study. 

C. Running in the U{2)l x U{2)r linear sigma model 

In component fields (19) can be written as 

1- ^2 ™',.2^1.;, ^2 "^a2 



^u{2) =il'2[i0 - g{^5 - il5a5)]ip2 + T^idf.'pf - -ifv'^ + T^id^a) 



a 



--y + a^f-^i^'a^-i^afl (20) 

where we introduced the following notations: (/?5 = CTq + TiUi, a^ = — ttq — TjTTi, with t^ being the Pauli matrices, 
A = Ai + A2/2, Lp = (cr, TTj), and a = {—ri,ai). The squared masses are m^^/^ = ^ T c, where m^ — M^ + Aij/^ and 
c=Cy{d.{5)). 

The running of the couplings are determined by the beta- functions given in Appendix A, under equation (A32). 
By solving (A32a) the running of g can be obtained explicitly: 

9 (P) = r-27-^ -2 = TT' (21) 

where Ag = /i§ exp [l27r^/(5.g^(/io))] • This has an UV Landau pole, while it goes to zero when fi ^ 0. 
For the other two equations, (A32b) and (A32c), we introduce the ratios: 

2A A2 

U = -^, U2 = ^J, 

9 9 

and a new function X{p) monotonous in ^ which satisfies 

1 dX 1 



(22) 



Using (21), this equation has the solution 



^(A*)-^ln 



In (Ao/^io) 



(23) 



ln:2i=in:^ e-lt^-^ -1 =ln^ 1-elt^-^ . (24) 



g'^ din 11 At:^ 



1^0 ,„ Ao ,^__iox i^_i ^0/-, mx 



ln(Ao/Ai) 

Here we have chosen the condition A(/io) = 0, where /iq is chosen to be the strange quark mass nis- Then we find: 

du 2 14 

— — = 4u + 3uu2 + 3m2 - 4 — —u, 
oX 3 

— =3mM2+M2 -4- — M2. (25) 

Conform subsection IV A, the phenomenologically motivated initial conditions are U2 = 2 and u = 2(^+1) = 4 ... 10, 
at the scale of the s-quark. The solution of (25) is depicted in Fig. 2. As this plot also demonstrates, for a wide 
range of /i, U2 stays in the interval [1,2], while u decreases continuously as we lower the scale fx from /ig = "^s down 
to /z < vris- Sooner or later (depending the initial conditions) u crosses zero which signals the instability of the theory. 
Probably it means that in order to maintain stability higher order corrections are needed. 
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FIG. 2. Running of the ratios u = 2\/g^ and U2 = ^ij 9^ starting from different initial conditions. 

D. Running in the SU{2)l x SU{2)r linear sigma model 

If we are well below the nia scale, then we can use the model containing the cr — tt sector of (20) and the u, d 
constituent quarks. This model is the SU{2)l x SU{2)ii linear sigma model defined by 



1 777 A 



(26) 



where ip2 and m^ were defined in the previous two subsections. The parameters of the model are defined at scale 
Too, so we have to apply renormalization group running to find the values of the coupling at the phase transition 
temperature T ^ m„. 

The RG equation are determined in Appendix B, eq. (B14). The running of g can be solved: 



g'ip^) 



g^M 



24^^ 



1 + ^i^ln^ 51n-^ 



(27) 



where Aq = /i§ exp[— 247r^/(55'^(/xo))]. This has an IR Landau pole, while it goes to zero when /i ^ 00. The definition 
of Ao is RG invariant. Comparing (27) with (21), and taking into account that the change of scaling is at jiq ~ nia 
we find 






(28) 



For the running of A we again introduce again X{p,) defined in (23) where now g is the coupling of the SU{2)l x 
SU{2)ji linear sigma model. Using (27) we obtain 



x{^l) 



■In 



ln(A^/Ao) 
ln(/io/Ao) 



In ^— = In -^ I 1 - e 3 

M ^0 



We introduce the same ratio as in (22): 2A = ?ig^, and find 



du 
dX 



3u^ + ?i - 4. 



The solution of this equation reads 



M-1 



Uq-1 



1 



ln(Mo/M) V^^^ 



3u + 4 3uo + 4 V ln(/^o/Ao) 
This equation has a fixed point at m = 1, as it is shown in Fig. 3 
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FIG. 3. Running in the SU{2)l x SU{2)ji linear sigma model showing the presence of a fix point a.t u = 1. 



In the SU{2)l X SU{2)ji linear sigma model we can follow the running of the light mass, using (B14). The same 
running is true for the tree-level sigma mass m'^ = — 2?Ti^. Then we find 



With the help of (27) and (31) we obtain 

mKfi) ^ m^(^o)exp 

21/5 



1 



127r2 



In- 



rf(ln(M'^o)).9'(A*')U^(A*')-l 



7g ifJ-a) /■ ™ ds l + 6Yoz{s) 



247r2 7o l + bg^fio)s 1 - 3Yoz{s) 



where z{s) = 1 + 



ln(^o/Ao) 



b = 5/(127r2) and Fq = (m(mo) - 1)/(3u(a*o) + 4). 



(32) 



(33) 



V. THERMODYNAMICS OF THE TRICRITICAL POINT 



The one-loop study of the tricritical point in the SU{2)l x SU{2)ii linear sigma model was done in [18] using an 
expansion in the number of flavors. We quote below the equations (12) and (14) of that work which determine the 
position of the tricritical point in the fiq — T plane. Using the present notation for the couplings these equations read: 



ml+T' 



, 2 , 3.g- 2 
-+g +^a 

3 TT^ 



-0, 



So* 



(34) 



where /? — 1/T is the inverse temperature and a — j3^q with fiq the quark baryon chemical potential. The function 
F reads 



d 



^(a) = 1 - 7£ + ln2 - — [Li,(-e") + Li,(-e-")] 



d. 



(35) 



This is a monotonously increasing function of its argument, T{Q) ~ 1.5675. 

We choose the scale fi — e^T, where ^ is a number of 0(1). Then the logarithm yields ^, which effectively modifies 
F ^^ F = F — i^. In case of spontaneous symmetry breaking m^ < 0, and it is useful to rescale all the masses with 
the tree level sigma mass rn^ := — 2m2 at scale ^. Then the equations to solve will be 



l = TV{f^) 



u(ri) _ 3 9 
^ + 1 + —a^ 

6 TT^ 



6 

(m) = ^5^(A*)-^(a)- 



(36) 
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FIG. 4. The dependence of Oc {fJ-q/T at the TCP) on the parameters for U2(ms) = 2. The left panel is obtained with 
u(ms) = 10. The value of Ao is fix along the lines of the surfaces directed towards the origin of the rUs — ma plane. For further 
information see the main text. 



In the complete problem therefore there are 5+2 parameters: at the UV scale tUs 
also we have /i and T at the IR scale. The light mass m^ or the corresponding jUa 
TCP a = Piiq and T can be determined as functions of the UV parameters: 

Go : ms,ma,g'^,u,U2 i-^ ac,Tc. 



we have nis, raa, g^ , u, U2, and 
is used as a mass unit. At the 



(37) 



The final output of the investigation should be, of course adms) and Tc{ms). But as we just have seen, even in 
the chiral limit of niud — the problem is five dimensional instead of one dimensional. For a sensible prediction we 
have to say something about the strange mass dependence of ma,g^,u,U2 ~ these functions should come from the 
underlying theory, now QCD. Since we do not have this information, we have to assume something sensible. 

In the light of the previous subsections we make some approximations: we can fix U2{ms) = 2 and for u we consider 
two cases: m(tos) — 4 and u{'ms) ~ 10. The remaining function 



g : ms,ma,g i-^ ac,Tc 



(38) 



can be plotted as shown in Fig. 4. The detailed numerical strategy to solve the system and obtain this plot is given in 
Appendix C. Fig. 4 shows surfaces in the nis, Jn^, g^ parameter space leading to some fixed value of a^. The piq/T — 
critical surface is a limiting one, in the sense that surfaces with [iqjT > all lie on one of its side, they never cross 
each other. Moreover the normal vector of the surface pointing to positive fXq/T always have negative Trig component 
- in this sense we can say that going on the direction of the largest fJ.q/T change, the surfaces bend downwards in mg. 
In order to show that different values of u do not change qualitatively the result we plot the Hq/T — 0.8 surface 
obtained using u{ms) = 10 and u{ms) — 6, and rescale the g^ values of the latter by a factor of 1.5. The two surfaces 
can be seen on the right panel of Fig. 4. As the plot shows, the surfaces have the same characteristics. 



A. Implication of the results for QCD 



In real QCD we cannot change ma, tUs and g^ independently. If we knew the nis dependence of iria and g^ then we 
would have a curve in the iris ~ raa ~ g^ space parametrized by mg . This line would go through the critical surfaces 
characterized by fix [iqjT , and then we could determine the ms{^q/T) function. Since we do not have any information 
on the rris dependence of the parameters we explore several possibilities by fixing the value of one of the parameters. 

For a constant value of g^ the Ws-dependence of the critical jiq/T is shown in Fig. 5. One can see that the behavior 
of this curve depends strongly on how rria depends on rris. Characterized by dma/dms at iiq = 0, there is a limiting 
value, and tricritical curves with smaller value of dma/dms bend downwards (negative curvature), for larger values 
they bend upwards (positive curvature). 

The standard characterization of the behavior of ms{fj.q) near /ig = is through the Taylor series [1]: 
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FIG. 5. Left panel: The ms-dependence of the critical iiqjT obtained for u = 10 and g^ = 3. The labels on the curves indicate 
the value of dma/{d,ms)(fi — 0). Right panel: The first two coefficients of the Taylor expansion in Eq. (39). 

The first two nontrivial terms C2 and C4 are shown in Fig. 5 in case of a constant g^. The singularity corresponds 
to that value of the mg for which the curvature changes sign. It is remarkable that by changing continuously from 
negative to positive curvatures (c2 values) we have to go through a singularity. 
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FIG. 6. Effect of a /:i-dependent ma on the m,s-dependence of the critical f^q/T obtained for u = 10. From the two pairs of lines 
having the same values of A/rricr and ma/rricr, the one at the right is obtained by solving (36) with a /i-dependent nia- 

To have a hint on which curve could be the physical one we recall that the anomaly is mainly a gauge effect, 
connected to the presence of the instantons [32, 33]. This suggests that The dependence of nia on ms should be quite 
small, so the physical line is near to dma/dnis ~ 0. Note that all tricritical surfaces in this model with nia =constant 
bend downwards. 

We can compare our results with the numerical findings of [1]. They found C2 = —3.3(3), C4 — —47(20) with Nf — 3 
degenerate quarks. The C2 curve on Fig. 5 reproduces this value at dnia/dms = —0.58 with a corresponding value of 
C4 = —33.0. This is in the same order of magnitude as in the MC simulation. 

In Ref. [25] it was shown that if the strength of the U{\) anomaly, parameter C in the Lagrangian (4), is made 
/ig-dependent then the critical surface can have a nonmonotonic shape. Since C influences ma, we can observe the 
same effect by considering the following dependence of rua on /i^, when solving (36): 



maifiq 



= 1 



ma{0) 



- 1 



2 / 2 



(40) 



For most of the curves with this chemical potential dependent anomaly are very similar to what is obtained earlier, 
however one can observe more exotic behavior, too, for some cases. For fiq.o/m^ = 0.17 the result is shown in Fig. 6: 
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here the surface starts to bend upwards, and later it turns back. This behavior can be very mild (as in our example 
with A/wg- = 208). In this case coarse lattice measurement would only detect the negative curvature, a high precision 
lattice measurement is necessary to reveal the positive curvature near fiq = 0. It is interesting that, considering the 
global behavior of the curve, the coarser lattice would give a more reliable result in this case. 

VI. CONCLUSIONS 

We discussed the behavior of the line of tricritical points (TCP) in the chiral limit {iriud — 0) of the U{3)l x U{3)]i 
quark model. We assumed that the value of the strange quark mass, where the TCP hits the fi — line is much larger 
than the critical temperature Tc- This is a good approximation in QCD, where the critical temperature is of order 
160 MeV, while the constituent strange quark mass is about 450 MeV already at the physical point, and we expect 
that the second order line reaches the chiral line {iriud = 0) at much higher rus masses. We also assumed that the 
r] — ao meson sector, which is heavier than the a — n sector because of the anomaly, is also much heavier than Tc- 
This is again plausible, since already at the physical point ma ~ 1 GeV. 

Under these circumstances the strange and the rj—aQ sector decouples from the point of view of the thermodynamics, 
which is completely determined by the light degrees of freedom, the a — tt sector. The only way how the heavy sector 
can influence the thermodynamics is through the values of the parameters of the effective theory. To achieve this goal 
we have to follow the running of the different parameters as well as the degrees of freedom from the heavy scales down 
to the thermodynamic scale. This can be performed by following the renormalization group (RG) flow with given 
degrees of freedom, and determine the parameters of the effective theories by matching when the degrees of freedom 
change. The former yield logarithmic dependence on the heavy scale, the latter effect is power-suppressed. Therefore 
in this work the RG flow is determined at one-loop level, and the matching is kept at tree level. 

With the decoupling of the strange and 77 — oq meson sector, respectively, there are two stages of effective models 
in the C/(3)l x C/(3)fl linear sigma model: the first is the U{2)l x C/(2)fl, while the second is the SU{2)l x SU{2)ji 
linear sigma model. We determined the corresponding beta functions in these models, and solved the RG flow down 
to the scale of the temperature T. The thermodynamics is determined at one-loop level in [18] - we now included the 
running coupling constant in the result. 

As a result we can determine the free energy for any given parameter sets, and we can determine those points where 
tricritical points (TCP) are located. With some plausible assumptions, the TCP's with fixed a = fJ-q/T (where /i^ is 
the quark chemical potential) form a surface in the m,, — nia — g^ space in the U{3)l x 1/(3)^ linear sigma model. 
Surfaces for different a never cross, therefore the a = surface is a limiting surface. 

If we want to draw consequences for QCD, we have to specify how nia and g^ depend on mg in the chiral u, d 
regime. Since it is not known, we explored several possibilities. Depending on the details, this model can describe an 
upward bending (positive curvature) surface or a downward bending (negative curvature) surface. Taking into account 
the explicit chemical potential dependence of the anomaly constant, the curvature we can change from positive to 
negative curvatures along the curve. In this case the value of the curvature ai fiq — would yield false information 
about the global behavior of the curve. 
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Appendix A: Renormalization group equations in the U{2)l x U{2)r model 

We start from the renormalized Lagrangian of (20), and introduce the countcrtcrm Lagrangian which in Fourier 
space reads 



^-l^ik'-Sml)^+^-^aik'-6ml)a-^-^{^' + ar-^-^'-'^' 



6Cui2) = 5Z^ij2H2 - Sgip2{^5 - il5a5)ip2 + —^{k^ - Sm^^)^ + —a{k^ - Smi)a -—{ip^ + a^ - -^ {^^a^ - {^aY), 

(Al) 

ip/a 



where we used the shorthand dm?' , — 6m? T 5c and the observation that the wave function renormalization for the 

in n. I 



If and a sector is the same. 

The goal is to determine the counterterms at one-loop level. To this end we work at zero temperature in the 
symmetric phase. The method is to determine the expectation value of some physical observables and require finiteness. 
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In order to simplify the treatment we introduce a background field for the a field through the shift cr — > a + x. The 
expansion in x is used for zero momentum external legs. The new Lagrangian is obtained from (20) reads 

l^vm =-^-' - 7-' - --K + ^-') + ^2(^^ - m^)^2 + \{d,^f + i(a,a)2 - ^a^ - ^nf ^rj' ^a' 



2 - 4~ -~v->p . ..^ , : -r.y^. ..-^,^. 2^ '^'^^ 2^ '^ ' 2 2 ^ 2 ' 2 

-g'ip2{'P5 - »75a5)V'2 - Aa;CT(v3^ + a^) - X2xaaf - X2xr]TTiai - jiv'^ + a^T ^ ['^'^a^ - {^o,f] , (A2) 

where 

m^ — gx, rn^ = m^ + 3Xx^, rn^ — m^ + Xx^ , m^ — m^ + Xx^, m\ = m^ + Xx^ + X2X^ , (^3) 

with m^ and m^ defined below (20). 

1. The fermionic wave function and g renormalization 

We calculate on the x background the fermion self-energy S^ = i (Ttp2'>p2) ■ Introducing the notation 
and using standard Feynman rules we find 



S^(fc) = -SZ^^ + xSg - ig^ f \{iG„{p - k) + 3iGa{p - k))ig{p) - {iGnip - k) + SiG^ip - k))-f5iGiph5\ , (A4) 
where G and Q are the bosonic and fermionic propagators defined as 

Gip) = - ^y--, Qip)= /^r . , (A5) 

with the corresponding masses. 

Since the integral has mass dimension, the part proportional to |^ or a; are dimensionless which means that they are 
at most logarithmically divergent. Therefore the masses should be taken into account only through a Taylor expansion. 
But, the expansion in the bosonic masses yields rn^ ^ x^ terms, which are convergent, so we can forget about the 
bosonic masses. The same is true for the fermionic mass in the denominator. What remains for the divergent piece is 



E^-(fc) = -5Z^f + x6g + 42.g2 / G{p - k) [g{p) - 75^(^)75] ■ (A6) 



i-p 

In the numerator we find ^ + ?7i — 75(^ + m)^s ~ 2^, which results in 

bg = 0. (A7) 

For 8Z^ we have to calculate the remaining integral. Doing this with standard techniques (cf. for example [34] ) using 
cut-off regularization we find 

and so, the counterterm ensuring the finiteness of E^(fc) is 

^^^--^2^^^- (A9) 

2. The bosonic wave function and A renormalization 

We calculate next the a self-energy on the given x background. We find 

S„(/c) = -6Zk^ + (5m2 + UXx^ + 3A / iG„{p) + 3A / iG^P) + A / iG^{p) 

J p Jp Jp 

+3(A + A2) f iGa{p) - tg" f Tr [g{p - k)g{p)] + ISiX'^x^ I G^p- k)GAp) 

Jp Jp Jp 

+6iAV I G^{p-k)G-,{p) + 2iX^x^ f G^ip-k)Gr,{p) + 6iiX + X2)^x^ j Ga{p- k)Ga{p). (AlO) 

Jp Jp Jp 



5Z = -^\n—, <5m2 =-_ln — [6Am2+(4A + 3A2)m2]. (A13) 
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The minus sign is because the fcrmionic bubble involves a closed fermion loop. To determine 5Z and 5m^p we need 
only the x = sector: 

S,(fc, a; = 0) = -SZk' + 8m\ + 6A / zG^(p) + (4A + 3A2) I iGM - *ff' / Tr {Q{v - ^)Q{v)\ ■ (All) 

Jp Jp Jp 

After evaluating the integrals, without writing the A^ corrections we find 

Sr (fc, . . 0) = -SZk^ + Sml + ^^mlln§ + 1^ ruling ^^k^ln^. (A12) 

Therefore 

9' A2 . 1 A2 

12'K-' II'' ^ Idtt^ 11 

For the determination of SX we need the self-energy at fc = 0. After evaluating the integrals we find 
J:f^k = 0)= 5ml - TTT^ In ^ [OAm^ + (4A + 3A2)m2] + 3S\x^ _ ^ In ^ [13A2 + 3(A + As)" - Ag*] . (A14) 
We obtain for Sml the previous result given in (A13), and we also have 

6X = -^ In ^ [13A2 + 3(A + A2)' - 4/1 . (A15) 

3. Renormalization of A2 

We apply the procedure above, but now for the a self-energy. Since the wave function renormalization is the same 
as for ip we need only the k = case. We find: 

P 

2 
p J p J p 

+4^(A + \2fx' f G„{p)Ga{p) + i\W I G^{p)G„(p). (A16) 

J p J p 

After evaluating the integrals we find for the divergent pieces 
Y.'^'^ik = 0)= 5ml + ('5A + 5\2)x^ - — ^ In ^ [GAto^ + (4A + 3A2)m2 1 - -^ In ^ [iGA^ + I8AA2 + SA^ - 129^] . 



T.a{k = 0) = (5m^ + ((5A + 5\2)x' + 5A / iGa{p) + (3A + 2A2) / iG^{p) 

J p J p 

+A j iG^ip) + (A + A2) / zG^b) - ig"" f Tr [g(rtg(rt] 



(A17) 



For the a-mass counterterm we find he following expression 



1 , A 



2 



S^l = Y^^^ — [<^^^l + i^^ + 3^2)ml], (A18) 

which is the same as the expression for the a mass, with the m^ ^-> nia interchange. This shows that the sum and 
the difference of the masses are renormalized multiplicatively: since "t,^ , — m? =p c, so we find 



2A - 3A2 , A2 ^2 2 lOA + 3A2 , A2 

— -^— In— , 5m ^m — — ^ — In — . 



'^c = c— — ^In— , 5m ^m ^^_, In—. (A19) 



From (A17) we can also read off the counterterm for A + A2 



1 A^ 

5\ + 5\2 = 7 In -^ [IGA^ + I8AA2 + 5X1 - iV] • (A20) 



Comparing it with the expression of 5X given in (A15) we find 



5X2 = -^ In ^ [A2 (6 A + A2 ) - 4g4] . (A2 1) 
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4. P functions 



The bare-field Lagrangian reads 



^u{2),o = i'mii^ - 9a{'P05~il5aa5)]ip02 + ^{d^(po)'^ + -^{d^aa) 



(A22) 
where aU the fields and couplings are bare. The bare couplings are RG invariant, since they depend only on the 
rcgularization: 



dgo 



dm^^ _ dml^ 



dXf) 



dX 



02 



dln/i dln/i rfln/i dln/i dln/i 



0, 



(A23) 



where 



d 



d\nfi d\nfi 



d d d 



'dg 



dX 



■d\2 



■Iv 






drn^, 



(A24) 

1/2 

To obtain the bare quantities from the counterterms, we first have to change to renormalized fields 1^02 — ZJ -02, 
(^0 = Z'^^'^V^ ^nd Oq = Z^/^a: 

Z Z 

£2 = ^2[Z^i$ - Z^Z'^/^ga{(p5 - ij5a5)]'4'2 + -^{d^j^ipf + —{d^af 



Zm, 



0^^2 



Z?7l 



0a„2 -^ '^0/_2 , „2n2 ^ Ao2r 2„2 



-(^^+«^)^ 



ip a — if a) 



(A25) 



2 " 2 4 '" ' 2 

Comparing it with the renormalized Lagrangian defined as the sum of (20) and (Al) wc find Z^ = 1 + SZ^ and 
Z — 1 + 6Z, so that the relations between the bare couplings and counterterms read 

Z^Z^/'^go=g + Sg, Zml^ = ml + Sml, Zml^ = ml + Sml, 

Z^\o^\ + SX, Z^Xo2^X2 + 6X2. (A26) 

These relations can be inverted, and at one-loop level we obtain: 

'1 
—( 

2 

Xa = X-2SZX + SX, Xa2 = X2 - 2SZX2 + SX2. (A27) 

Perturbative hierarchy requires that when there is a In fi dependence in the quantity, then only the d/{d In fi) derivative 
acts on it. Then, using (A24) we find 



go= g- {-SZ + SZ,p] g + Sg, ml^ ^ ml - SZm^ + Sm^, ml^ ^ ml - SZml + Sml, 



dgo 
dln/x 



/?.- 



d 



and in consequence 



In a similar way we find 



^ dlnij 
d 



-SZ + SZ^ ) g - Sg 



= 0, 



dhifi 



-SZ + SZ^ ] g~Sg 



(A28) 



(A29) 



5 



7v = 



[SZm — 5m 



d 



la = 



91n/i 



[SZml - 5ml] , 



din ji 

I3x - TTT— [2<5ZA - ,5A] , P2 = ttt— [25ZX2 - 8X2] . 
a m yu a In ^ 

Using the expression of the counterterms determined in previous subsections of this sections we have: 



(A30) 



dSZ, 



f_ d6Z _ g^ 

din ji 27r^ ' 9 In /i Gtt^ 



V- 



dSg 
din 11 



= 0, 



d6X 
dlnfi 



1 



87r2 
9 In /i Stt^ 



13A2 + 3(A + X2f - 4/1 , 



86X2 
91n/i 



47r2 



[A2(6A + A2)-454 



9Sml/^ _ 1 



6Am2/„ + (4A + 3A2)m2/^ 



(A31) 
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With these expressions, we obtain from (A29) and (A30) the fohowing one-loop /3-functions: 

dg 5g^ 



dlnfi ^^ 127r2 
dA _ 1 



dln/i 
dX2 



13A2 + 3(A + A2)2 - 4/ - Is'A 



1 



dm? I 1 

l^ I a 



A2(6A + A2)-V- ^2^2 



d\ii\i 



87r2 



(6A-^5')rn^/„ + (4A + 3A2K/^ 



(A32a) 
(A32b) 
(A32c) 
(A32d) 



Appendix B: Renormalization group in the SU{2)l x SU{2)r linear sigma model 

We start from the renormaUzcd Lagrangian (26) and add to it the the foUowing countcrtcrm Lagrangian: 



Z , 



6m , 



6^SU(2) = V'[^i/'*^ ^ ^9{<^ + JTiTTi)] V' + —{d^ipf -^ip- 



5X 



(^2 j_ Jls^ 



(Bl) 



The goal is to determine the counterterms at one-loop level. To do this, we will follow the same strategy as in 
Appendix A. To facilitate the discussion we introduce again the background field x. After the shift a ^ o -\- x the 
Lagrangian reads 



£ 



SIJ{2) 



where 



V J2 



A 



1 



1 



-X - -X - ax{m^ + Xx ) + ■)p{i0 - m^)-)p - g-)p{a + iTi7ri)-)p + -{df_,a) ^a + -(9^7r) 

-Xxa{a^+7rf)-^ia^+nfr, 



m,^ = gx, m^ = m'^ + 3Xx'^ , m^ = m^ + Aa;^. 



— n 
2 

(B2) 
(B3) 



1. The fermion Avave function and g renormalization 

We calculate the fermion self-energy on the background x: 

Y,^{k,x) = -SZ^jt + Sgx - ig"^ / [iGaip - k)ig{p) - SiG^rip - fc)75«^(p)75] , (B4) 

Jp 

where Go-, G^, and Q are the bosonic and fermion propagators introduced in (A5) with the corresponding masses 
given in (B3). Taylor expanding the bosonic propagators in x and using that 75 anticommutes with all the other 
gamma matrices, we find 



S^(/c) ~ Sgx — SZ^jt + 2ig 



2^)- m^ 



Ip {{P - kY - m2 + ie)(p2 - m| + ie) ' 
After evaluating the integrals, finiteness of the result requires the following expressions for the counterterms 



A2 ^ 53 . A2 






(B5) 



(B6) 



2. The (J mass, wave function and A renormalization 



Next, we calculate the a self-energy: 



E^(fc) =-5Zk^ + 5ml + 3'5Aa^^ + 3A / iGAp) + 3A / iG^(p) - ig'^ / Tr [g{p - k)g{p) 

J p J p 

-MSaV f Ga{p-k)Ga{p) + 6iX^x^ j G^{p- k)G^{p). 

J p J p 



(B7) 
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The minus sign is because the fermionic bubble. After evaluating the integrals we find, neglecting the A^ corrections 

3A 



Ediv(^) = -dZk^ + 5ml + 3(5Ax2 - — ^ 

^ lOTT-^ 



^o 2 .X 2N, A' 9' fk^ o 2 2\, A2 3A2, A^ 
(2m2 +4Ax2)ln- + ^ (^ y + 3^^ ) In ^ - ^ In ■ 



^2 27r2 ^2 



Therefore, the expression of the counterterms is 

„2 ^2 



6Z = 



In- 



3A 



247r2 ^2 ' 



Sm,„ ~ -—^ m,„ In 



A2 



'^ 87r2 "^ ^2 



,,.2^1^,„A^ 



47r2 



2 ■ 



(B8) 



(B9) 



3. /^-functions 



We again use the fact that the bare couplings are renormalization group invariant: 

dgo drriQ dAg 



dlnfi dlnfi d In /i 



= 0, 



where 



This leads to: 



d 



dln/i 91n/i 



/3,^+/3. 



d_ 



-Iv 



dm^ 



91n/i 



SZ + -SZ^ ] g - Sg 



din 11 



[SZm, 



y Sm^l 



din 11 



[2SZX - SX] . 



Using the counterterms determined in the previous two subsections one has: 



ddZ 



V _ 5 



dSZ 



dSg 



9 



d6X -3A2 + c,4 



9 In /It 47r2 ' dln/j, 127r2 ' 91n/i 47r2 ' dln/j, 

Then, we find the following one-loop /^-functions 

dg ^ 5.g3 dX ^ 9X^ - 3g'^ - Xg^ 



27r2 



(91n/i 



47r2 



dliifi 



P9 



247r2' 



dln/i 



/3a 



67r2 



1 dm^, 
m^ d In ^ 



9A-g2 



(BIO) 



(Bll) 



(B12) 



(B13) 



(B14) 



Appendix C: Numerical strategy to solve the TCP equations 

For numerical purposes it is advantageous to choose a, g^ = g'^ip) and nia as parametrization variables. Then we 
can proceed as follows. From (36) we find 



6J"(a) 2 

.77-Z 



rji'2 __ 





1 


2g^ 


"S^ + 1 + ^a^' 



(CI) 



then from (27) 



Ao = A^^e '5s^(M) 





r,i- 247r^ 


2gHp) 


[»'(^?(")+l + 4a2] 



(C2) 



Once we know g^{n), T, Aq and u{ii) we can compute mo- by solving m^ = m^d^ = nia) equation using (33). Since 
now ma-{i^) = 1 is the mass scale, in view of (33) we have to solve 



rria- — cxp 



730^ ^'"^ 



ds 1 + GYzis) 

24tt^Jo l + bg'^{fi)s l-3Yz{s) 



(C3) 
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where Y = (u(/z) — l)/(3u(/z) + 4). If we know nic then wc can have g'^{raa) from (27) which can be kept fixed. 
From the running of g (27) and u (31) we find 



g^ima) 



2u{ma) - 1 _ 2w(/i) - 1 A , \n{ma/li)\^^'^ 



5^2 m„ ' 3w(ma) + 4 3m(m) + 4 V" ' ln(A*/Ao) 

127r'^ ^ 



(C4) 



From (28) wc find: 



^t = 2mlg\^,) 



g^{y)F{a) 



l + ^a' 



g 53^(m) , 



(C5) 



Having Ag, u{ma) and u{ms) we can use the solution of the U{2)l x U{2)fj RG running (25) to find Trig. Finally from 
TOs and g(ma) we compute from (21): 



/("Is) = 



g^{ma) 



1 -- — In — 



(C6) 
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